Theory of Electro-Optical Properties of Graphene Nanoribbons 
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We present calculations of the optical absorption and electro-absorption spectra of graphene 
nanoribbons (GNRs) using a n— electron approach, incorporating long-range Coulomb interactions 
within the Pariser-Parr-Pople (PPP) model Hamiltonian. The approach is carefully bench marked 
by computing quantities such as the band structure, electric-field driven half metallicity, and linear 
optical absorption spectra of GNRs of various types, and the results are in good agreement with 
those obtained using ab initio calculations. Our predictions on the linear absorption spectra for 
the transversely polarized photons provide a means to characterize GNRs by optical probes. We 
also compute the electro-absorption spectra of the zigzag GNRs, and argue that it can be used to 
determine, whether or not, they have a magnetic ground state, thereby allowing the edge magnetism 
to be probed through non-magnetic experiments. 

PACS numbers: 78.20.Bh, 78.67.Wj, 73.22.Pr , 78.40.Ri 
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I. INTRODUCTION 

Discovery of graphene 1 has stimulated intense re- 
search in the field from the point-of-view of both fun- 
damental physics, and promising applications^—. Of 
particular interest are recently synthesized quasi-one- 
dimensional (ID) nanostructures of graphene called 
graphene nanoribbons (GNRs) which have technologi- 
cally promising electronic and optical properties because 
of the confinement of electrons owing to the reduced di- 
mensions. As a result, numerous theoretical studies of 
electronic, transport, and optical properties of GNRs of 
various type have been performed over the years£~— . The 
structural anisotropy of GNRs must exhibit itself in an 
anisotropic optical response with respect to the photons 
polarized along the length of the ribbons [x polarized, or 
longitudinally polarized) as against those polarized per- 
pendicular to it (y polarized or transversely polarized), 
with GNRs being in the xy-plane. Despite its obvious 
importance, anisotropy in the optical response of GNRs 
has not been studied in any of the reported optical ab- 
sorption calculations, which concentrate only on the lon- 
gitudinal component of the spectra*^— . In this work 
we study this anisotropy in detail, and make predictions 
which can be tested in optical experiments on oriented 
samples of GNRs, and can serve as a means for their 
optical characterization. 

Electro-absorption (EA) spectroscopy, which consists 
of measuring optical absorption in the presence of a static 
external electric (E) field, has been used extensively to 
probe the electronic structure and optical properties of 
conjugated polymers and other materials^!. GNRs, being 
7r-conjugated systems, will also be amenable to similar 
EA probes, and, therefore, we have calculated the EA 
spectrum of zigzag GNRs (ZGNRs) in this work. ZGNRs 
have been predicted to possess a magnetic ground state, 
with oppositely oriented spins localized on the opposite 
zigzag edges of the ribbons^. Our calculated EA spectra 
of ZGNRs depends strongly on whether, or not, they 
exhibit edge magnetism, thereby, allowing its detection 



by optical means. 

Most of the theoretical approaches used to study the 
electronic structure of GNRs are broadly based upon: 
(a) tight-binding method^^ (b) Dirac equation ap- 
proach, derived using the linearity of the band struc- 
ture in the region of interest ^ (c) ab initio DFT 
and GW based approaches) 10 ! 11 ' 13 ! 15 and (d) Hubbard 
model based approaches But, it is obvious from 
the chemical-structure of graphene and GNRs that the 
electrons close to the chemical potential are itinerant ir- 
electrons which determine their low-energy excitations. 
In 7r-electron systems such as various aromatic molecules 
and conjugated polymers, it is well-known that role 
of electron-electron (e-e) interactions cannot be ignored 
when describing their electronic properties.- 23 Therefore, 
it is inconceivable that the long-range e-e interactions will 
be insignificant in graphene and related structures. The 
effective 7r-electron approaches such as the Pariser-Parr- 
Pople (PPP) model Hamiltonian^ which incorporate 
long-range e-e interaction, have been used with consid- 
erable success in describing the physics of 7r-conjugated 
molecules and polymers^ Computationally speaking, 
PPP model has the advantage of including the long-range 
Coulomb interactions of it electrons within a minimal ba- 
sis, thereby allowing calculations on such systems with 
limited computer resources, as compared to the ab initio 
approaches. Indeed, in our earlier works we have used the 
PPP model to extensively to study the electronic struc- 
ture and optical properties of finite 7r-electron systems 
such as conjugated molecules and oligomers at various 
levels of theory^ 5 - Therefore, in this work, we have de- 
cided to extend our PPP model based approach to study 
the physics of GNRs in the bulk limit. Because, to the 
best of our knowledge, this is the first application of 
the PPP model to the GNR physics, we have carefully 
bench marked it for quantities such as the band structure, 
electric-field driven half metallicity, and linear optical ab- 
sorption spectra against the published ab initio works on 
GNRs, and the results are in very good agreement with 
each other. 

The remainder of this paper is organized as follows. 



In the next section, we outline the theoretical aspects 
of our work. In section IIII1 we present and discuss our 
results. Finally, in section HVl we present our conclusions 
and discuss the directions for the future work. 



II. THEORETICAL DETAILS 

The PPP model Hamiltonian^ with one 7r-electron 
per carbon atom (half- filled case), is given by 

H = - %(4rC?V + c) a C ia ) + 

u^mtnn + ]T Vijim - - l) (l) 

% i<j 

above cj CT creates an electron of spin a on the p z orbital 
of carbon atom i, — c ia .Ci a is the number of elec- 
trons with spin cr, and rij = ^ CT rii a is the total num- 
ber of electrons on atom i. The parameters U and Vij 
are the on-site and long-range Coulomb interactions, re- 
spectively, while tij is the one-electron hopping matrix 
element which, if needed, can be restricted to nearest- 
neighbors (NN). On setting Vij — 0, the Hamiltonian re- 
duces to the Hubbard model. The parametrization of the 
Coulomb interactions is Ohno like 2 ^, 

Vi,j = u/K id {i + omnRlj) 1 / 2 , (2) 

where, Ki t j depicts the dielectric constant of the sys- 
tem which can simulate the effects of screening, and Rij 
is the distance in A between the i-th and the j-th car- 
bon atoms. The Hartree-Fock (HF) theory for periodic 
one-dimensional systems, within the linear combination 
of atomic orbitals (LCAO) approach is fairly standard, 
and we have implemented both its restricted (RHF) and 
unrestricted (UHF) variants. The lattice sums are per- 
formed in the real space by including a large number 
of unit cells, and integration along the Brillouin Zone 
(BZ) was performed using the Gauss-Legendre quadra- 
ture approach 2 —. The convergence with respect to the 
numbers of unit cells included in the lattice sums, as 
well as fc-points used for BZ integration, was carefully 
checked. 

Our calculations, to the best of our knowledge, are the 
first applications of the PPP model to GNRs in the bulk 
limit; therefore, it is important to obtain a suitable set of 
Coulomb parameters for these systems. In our previous 
calculations on conjugated molecules and polymers^, we 
used two sets of Coulomb parameters namely (a) "stan- 
dard parameters" with U = 11.13 eV and Kij = 1.0, 
and (b) "screened parameters" with U — 8.0 eV and 
Kij = 2.0 (i ^ j) and K^i = 1, proposed initially by 
Chandross and Mazumdar to study phenyl-based con- 
jugated polymers^ In the absence of extensive exper- 
imental data, we adopted the criterion of good agree- 
ment between the ab initio GW band gaps of armchair 
GNRs (AGNRs)^ 3 . and our PPP band gaps, to choose 
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Figure 1: The structures of (a) a ZGNR and (b) an AGNR. 
The ribbons are assumed to lie in the xy plane, with the 
periodicity in the x direction. 



the Coulomb parameters. The tuning of the parameters 
was done for AGNR- 12 (AGNR- AT^ , denoting an AGNR 
with Na dimer lines across the width), and with a mod- 
ified set of screened parameters (U = 6.0 eV,K;j = 2.0 
(i j) and Ku — 1), and NN hopping matrix element 
<i = —2.7 eV. As a result, good agreement was obtained 
for AGNR- 12 between the PPP band gap (1.75 eV) and 
the corresponding GW value of Yang et a/. 13 . Therefore, 
we have decided to use these modified Coulomb param- 
eters throughout these calculations, with the aim that 
they will incorporate the GW-level electron-correlation 
results implicitly in our results. 



III. RESULTS AND DISCUSSION 

The schematic structures of AGNRs and ZGNRs stud- 
ied in this work are presented in Fig. [TJ Next, we present 
the results of our PPP model based calculations on vari- 
ous quantities, for AGNRs and ZGNRs. 



A. Band Structure 

In Fig. [5] (a) we present the band structures of AGNR- 
11 obtained using the Hubbard model with U — 2.0, and 
the PPP model. At the tight-binding level all the AGNRs 
with Na = 3p + 2 (pa positive integer) are predicted to 
be gapless. However, ab initio DFT calculations predict 
all types of AGNRs to be gapped, including those with 
N A = ip + 3ili 2 .. Our RHF calculations are m agree- 
ment with the DFT results, and also predict all families 
of AGNRs, including A^4 = 3p + 2 to be gapped, as is 
obvious from our PPP results for AGNR- 11 presented in 
Fig. [5] (a). The noteworthy point is that the Hubbard 
model, with the currently accepted values of U predicts 
a negligible gap for Na — 11 [cf. Fig. [5]), a result in 
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Figure 2: (Color online) Band structure near the Fermi en- 
ergy (Ef = 0) of: (a) AGNR-11 obtained using the Hubbard 
Model (black solid line), with U = 2.0, and the PPP-RHF 
approach (red broken line), (b) ZGNR-12 obtained using the 
PPP-RHF approach for the non magnetic state (black solid 
line), and the PPP-UHF approach for the magnetic state (red 
broken line) in which the bands of up and down spins are de- 
generate, (c) ZGNR-16, obtained using PPP-UHF model, in 
the presence of a lateral electric field of 0.1 V/A so that the 
up-down degeneracy is lifted (red broken/black solid bands 
represents up/down spins) with = o.ll eV (magnified 

in the inset), and E f iown) = 0.97 eV. 



complete disagreement with the DFT, and our PPP re- 
sults. Thus, from this case it is obvious that for AGNRs, 
long-range Coulomb interactions as included in the PPP 
model play a very important role of opening up the gap 
for the Na — 3p + 2 case. Our PPP value of the band 
gap 1.06 eV of this AGNR is again in excellent agreement 
with the ab initio GW result reported by Yang et alr^. 

The case of the ground state of ZGNRs is an in- 
teresting one with several authors reporting the exis- 
tence of a magnetic ground state, with oppositely ori- 
ented spins localized on the opposite zigzag edges of the 
ribbons^, a result verified also in several first principles 
DFT calculation o 11 ' 15 . We investigated this in our PPP 
model calculations by using the RHF method for the non- 
magnetic state and the UHF method for the magnetic 
one, and the results are summarized in Table |H We find 
that for a ZGNR of width Nz (Nz = number of zigzag 
lines across the width), ZGNR-iVz in short, the total 



Width 


Total Energy (eV) 


Band gap (eV) 


N z 






AE 




/rim 


4 


-23.059 


-23.261 


-0.025 


0.524 


2.414 


6 


-35.559 


-35.825 


-0.022 


0.336 


2.005 


8 


-48.103 


-48.403 


-0.019 


0.246 


1.694 


12 


-73.237 


-73.570 


-0.014 


0.161 


1.287 


16 


-98.417 


-98.743 


-0.010 


0.046 


1.037 



Table I: Variation of total energy/cell and the band gaps of 
ZGNR with the width of the ribbon, computed using the PPP 
model. 



energy/cell of the magnetic state (E m ) is lower as com- 
pared to that of the non-magnetic (E nm ) one, and the en- 
ergy difference/atom between the non-magnetic and the 
magnetic states (AE — (E m — E nm )/N at , N at =number 
of atoms in the unit cell) decreases with the increasing 
ribbon-width, consistent with the non-magnetic ground 
state of graphene. The band gap for the magnetic state 
(Eg 1 ) is much larger than that of the non- magnetic one 
(Eg 1 " 1 ). The non-zero gaps obtained for the non-magnetic 
states of ZGNRs is an artifact of the RHF approach. The 
band structures of the magnetic and non-magnetic states 
of ZGNR-12 computed using the PPP model are pre- 
sented in Fig. [2] (b) , and it is obvious that, for the 
magnetic case, our results are qualitatively very similar 
to the reported ab initio band structure a 11 ' 15 . Quanti- 
tatively speaking, for ZGNR-8, we obtain E™ = 1.70 
eV, which is higher than the reported GW value of 1.10 
eV— . Our band gap for AGNR-11 was in excellent agree- 
ment with the GW value, but that is not the case with 
ZGNRs. We believe that it could possibly be because: (a) 
our Coulomb parametrization was based upon ab initio 
GW results^ on an AGNR, and (b) electron-correlation 
effects are stronger in ZGNRs as compared to AGNRs, 
and the HF approach adopted here ignores those effects. 

In a pioneering work Son et a/.—, based upon ab initio 
DFT calculations, predicted that in the presence of a lat- 
eral electric field, ZGNRs exhibit half-metallic behavior 
leading to their possible use in spintronics. They demon- 
strated that for the field strength 0.1 V/A, the gap for 
one of the spins of ZGNR-16 will close, leading to metal- 
lic behavior for that spin orientation. In Fig. [5Jc) we 
present the band structure of the same ZGNR exposed 
to the identical field strength, calculated using the PPP 
model, and the tendency towards half-metallicity is ob- 
vious. While the band gap in the absence of the field 
was 1.037 eV, in the presence of the field up-spin band 
gap is reduced to 0.11 eV, while the down-spin gap de- 
creases to 0.97 eV. Therefore, considering the fact that 
our PPP model based approach does not incorporate 
electron-correlation effects, its quantitative predictions 
are in very good agreement with the ab initio ones 10 , 
and thus it is able to capture the essential physics of the 
electric-field driven half-metallicity in ZGNRs. 



B. Optical Absorption 

Next we present the linear optical absorption spectra of 
GNRs, computed within the PPP model. The optical ab- 
sorption spectrum for the x-polarized (y-polarized) pho- 
tons is computed in the form of the corresponding com- 
ponents of the imaginary part of the dielectric constant 



tensor, i.e. 



,(2) 



vvM). 



7T / a 



using the standard formula 

|< c (A0biKfc)>l 2 



-n/a {(Ecv(k) 



Ah 

■M 2 + 7 2 }^(fc) ' 

where a is the lattice constant, pi denotes the momen- 
tum operator in the i-th Cartesian direction, u> repre- 
sents the angular frequency of the incident radiation, 
E C v(k) = e c (k) — e v (k), with e c (k) (e v (k)) being the con- 
duction band (valence band) eigenvalues of the Fock ma- 
trix, 7 is the line width, while C includes rest of the con- 
stants. We have set C = 1 in all the cases to obtain the 
absorption spectra in arbitrary units. The components of 
the momentum matrix elements (c(k)|p|w(k)) needed to 

(21 

compute e H (w), for a general three-dimensional system, 
can be calculated using the formula,— 



(c(k)|pKk)) = ^< c (k)|v k ff(k)Kk)) 



im (e c (k) - e„(k)) 



c(k)|dKk))(4) 



where mo is the free-electron mass, Vk-ff(k) repre- 
sents the gradient of the Hamiltonian (Fock matrix, in 
the present case) in the k space, (e(k)|d|u(k)) denotes 
the matrix elements of the position operator d defined 
with respect to the reference unit cell, and accounts 
for the so-called intra-atomic contribution. 29 Note that 
Eq. 2] can also be used to compute the matrix element 
(c(k)\py\v(k)) needed to compute the absorption spec- 
trum for the y-polarized light for GNRs (which are peri- 
odic only in the x direction) , by setting the first term on 
its right hand side to zero, because for a one-dimensional 
system periodic along the x direction, the Hamiltonian 
has no k y dependence. (c(k)|Vk-ff(k)|u(k)) for the case 
of GNRs is obtained easily by calculating the numerical 
derivative of the Fock matrix at various fc-points of the 
one-dimensional Brillouin zone. For the d operator, the 
usual diagonal representation was employed. The cal- 
culation of the absorption spectra of the GNRs for the 

(2) 

y— polarized photons (ey y J (uj)), to the best of our knowl- 
edge, has not been done earlier. Because such transverse 
excitations do not couple to the photons polarized along 
the x direction, they have also been called "dark excitons" 
in the literatur e) 13 ! 15 . 

The optical absorption in AGNRs has been studied ex- 
tensively by ab initio approaches in recent work o 12 i 13 ' 16 . 
In Fig. Ela) we present the optical absorption spectrum 
of the AGNR-11. If £ mra denotes a peak in the spectrum 
due to a transition from m-th valence band (counted from 
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Figure 3: (Color online) Imaginary parts of the dielectric con- 
stant (eS(tu) in black solid, and e$y(u)) in red broken lines) 
computed using the PPP model, and modified screened pa- 
rameters, for: (a) AGNR-11, (b) ZGNR-8, with a magnetic 
ground state. Labels of the peaks denote the bands involved 
in the transition (see text for an explanation), and a line width 
of 0.05 eV was assumed throughout. 



top) to the n-th conduction band (counted from bottom) , 
the peak of e ( £}(uj) at l.leVis S u , at 3.1 eV is £ 22 , at 3.8 

eV is £ 33 , and at 5.8 eV is S 44 . The peaks of e^fw) at 2.1 
eV and 5.6 eV both correspond to S 12 and E 21 . The re- 
markable feature of the presented spectrum is that owing 
to the symmetry of the AGNRs, the peaks correspond- 
ing to x- and y-polarized photons are well separated in 
energy, and their relative intensities can be measured by 
performing experiments on oriented samples. On com- 
paring our PPP spectrum (e^x (^)) with the ab initio 
GW spectrum of Yang et al.—, we note that the locations 
of the first peaks close to 1.1 eV are in excellent agree- 
ment with each other. However, our calculations predict 
several higher energy peaks with significant intensities 
located around 3.0 eV absent in the GW work. Further- 
more, we also predict the intensities of the y-polarized 
peaks, which was absent in the work of Yang et ali^. 
In Fig. EJb) we present our calculated optical ab- 

(21 (2) 

sorption spectrum (e xx '(cj) and e yy (w)) for the ZGNR-8. 

The peaks in e x 2 j(w) are located at 1.7 eV (S 11 ), 2.9 eV 
(S 12 + E 21 ), at 4.0 eV (S 22 ), while the prominent peaks of 

4 2) H are at 1.7 eV (S 11 ) and 2.9 eV (E 12 + S 21 ). The 
noteworthy point is that most of the prominent peaks 
have mixed polarization characteristics, unlike the case 
of AGNRs. This is because of the fact that for mag- 
netic ground states, the reflection about the rrz-plane is 
broken, leading to mixed polarizations. This is an impor- 
tant result which can also be tested in oriented samples 
of ZGNRs. Our PPP optical absorption spectrum of this 



ZGNR compares qualitatively well to the GW spectrum 
computed by Yang et al. 15 , although our peaks are con- 
sistently blue-shifted compared to the GW result, due to 
the corresponding disagreement in the band structure. 
Moreover, Yang et alJ£ did not compute the peak inten- 
sities for the y-polarized photons. 

C. Electro- Absorption 

In Figs. 2] we present the EA spectrum of ZGNR-8 
computed as the difference of the linear absorption spec- 
tra with and without an external static E-field of strength 
0.1 V/A along the j/-axis. In Fig. 0Ja) we present the 
EA spectrum for the non-magnetic ground of ZGNR-8, 
computed using the PPP-RHF approach. Without the 
external E-field, the En transition is disallowed for the 
non-magnetic state of such a ZGNR for the x-polarized 
light due to symmetry selection rules 14 . However, in the 
presence of the field, due to the broken symmetry, this 
transition becomes strongly allowed leading to a very 
strong peak in the EA spectrum. Fig. HJb) portrays the 
EA spectrum of the same ZGNR for the magnetic ground 
state, and, here the physics of half-metallicity manifests 
itself in that one observes two energetically split peaks 
corresponding to two different Entransitions among up- 
and down-spin electrons. Thus, our calculations predict 
that the EA signal is different for the ZGNRs depend- 
ing on whether they have a magnetic or a non-magnetic 
ground state, a result which can be used to determine 
the nature of the ground state of ZGNRs using EA spec- 
troscopy. 



It will also be of interest to perform similar studies on 
bilayer and other multilayer GNRs, to investigate how 
various properties of the ribbons evolve, as the number 
of layers are increased. Of particular interest is the case 
of multila; " T ~ " e nature 

of edge m; e J ^ will 
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Figure 4: (Color online) Linear absorption spectrum (black 
solid) and electro absorption (red broken) of ZGNR-8 for pho- 
tons polarized along the x axis for: (a) non magnetic ground 
state, and (b) magnetic ground state. A line width of 0.05 
eV was assumed throughout, and the bands involved in the 
electro-absorption peaks are indicated. 



IV. SUMMARY AND OUTLOOK 

In summary, we have used a PPP model based 7r- 
electron approach, incorporating long-range Coulomb in- 
teractions, to study the electronic structure and optical 
properties of GNRs in the bulk limit. In particular, we 
computed the optical absorption spectra of GNRs for 
transversely polarized photons, in addition to the lon- 
gitudinal ones, thereby allowing us to investigate the 
anisotropic optical response of these materials. Our pre- 
dictions that for AGNRs longitudinal and transverse po- 
larized components will be well separated energetically, 
while ZGNRs will exhibit absorption with mixed polar- 
ization, can be tested in experiments on oriented samples. 
Furthermore, we also presented first calculations of the 
EA spectra of ZGNRs, and our results suggest a possibil- 
ity of an optical determination of whether, or not, they 
possess a ground state with edge magnetism. 



also be of interest to include excitonic effects in the op- 
tical absorption spectrum of ZGNRs so as to perform a 
complete comparison with the future experimental work 
on these systems. For the purpose, it is important to go 
beyond the Hartree-Fock approach and include electron- 
correlation effects. Work along all these directions is in 
progress in our group, and the results will be communi- 
cated in the future publications. 
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